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The heat kernel expansion for field theory at finite temperature is constructed. It is based on the 
imaginary time formalism and applies to generic Klein-Gordon operators in flat space-time. Full 
gauge invariance is manifest at each order of the expansion and the Polyakov loop plays an important 
role at any temperature. The expansion is explicitly worked out up to operators of dimension six 
included. The method is then applied to compute the one loop effective action of QCD at finite 
temperature with massless quarks. The calculation is carried out within the background field method 
in the MS scheme up to dimension six operators. Further, the action of the dimensionally reduced 
CO ' effective theory at high temperature is also computed to the same order. Existing calculations are 

reproduced and new results are obtained in the quark sector for which only partial results existed 
. up to dimension six. 
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I. INTRODUCTION 



The extension of field theory from zero to finite temperature and density is a natural step undertaken quite early 
QH>H>Q>IEI9' The interest is both at a purely theoretical level and in the study of concrete physical theories. At the 
■ theoretical level one needs appropriate formulations of the thermal problem, for which there are several formalisms 
^—1 \ available 0, as well as mathematical tools to carry out the calculations. From the point of view of concrete theories 
CO ■ a central point is the study of the different phases of the model and the nature of the phase transitions. That study 
applies not only to condensed matter theories but also to fundamental ones, such as the electroweak phase transition, 
of direct interest in early cosmology and baryogenesis 0, and quantum chromodynamics which displays a variety of 
Q H ' phases in addition to the hadronic one [3, [T3, HH LL3 • Such new phases can presumably be probed at the laboratory 
in existing (RHIC) 0] an d future (ALICE) facilities. Obviously one expects all these features of QCD at finite 
temperature to be fully consistent with manifest gauge invariance. As is well-known Lorentz invariance is manifestly 
broken due to the privileged choice of the reference frame at rest with the heat bath, however, gauge invariance 
• • ■ remains an exact symmetry. At zero temperature preservation of gauge invariance involves mixing of finite orders in 
perturbation theory. As will become clear below, compliance with gauge invariance requires mixing of infinite orders 
^ . in perturbation theory at finite temperature. 

The purpose of the present work is twofold. The first part (Section[n} is devoted to introduce a systematic expansion 
for the one-loop effective action of generic gauge theories at finite temperature in such a way that gauge invariance is 
manifest at each order. In the second part this technique is applied to QCD in the high temperature regime, firstly 
to compute its one-loop gluon and quark effective action (Section lllljl . and then to derive the Lagrangian of the 
dimensionally reduced effective theory (Section llV|) . Further applications can and will be considered in other cases of 
interest [l4| . 

The effective action, an extension to quantum field theory of the thermodynamical potentials of statistical mechanics, 
plays a prominent theoretical role, being directly related to quantities of physical interest. To one loop it takes the 
form cTrlog(.K'), where K is the differential operator controlling the quadratic quantum fluctuations above a classical 
background. Unfortunately, this quantity is afflicted by mathematical pathologies, such as ultraviolet divergences or 
many- valuation (particularly in the fermionic case). For this reason, it has proved useful to express the effective action 
in terms of the diagonal matrix elements of the heat kernel (or simply the heat kernel, from now on) (x\e~ rK \x), 
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by means of a proper time representation (see e.g. eq. (|2.17|l below) [110. Unlike the one-loop effective action, 
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the heat kernel is one-valued and ultraviolet finite for any positive proper time r (we assume that the real part of 
K is positive). A further simplifying property is that, after computing the loop momentum integration implied by 
taking the diagonal matrix element, the result is independent of the space-time dimension, apart from a geometrical 
factor. In practice the computation of the heat kernel is through the so called heat kernel expansion. This is an 
expansion which classifies the various contributions by their mass scale dimension, as carried by the background fields 
and their derivatives. This is equivalent to an expansion in the powers of the proper time r. In this way the heat 
kernel is written as a sum of all local operators allowed by the symmetries with certain numerical coefficients known 
as Seeley-DeWitt or heat kernel coefficients. The perturbative and the derivative expansions are two resummations 
of the heat kernel expansion. This expansion has been computed to high orders in flat and curved space-time in 
manifolds with or without boundary and in the presence of non-Abelian background fields 0, 0, 0, |2(J, 0, 152, • 

In order to apply the heat kernel technique to the computation of the effective action at finite temperature it 
is necessary to extend the heat kernel expansion to the thermal case. This can be done within the imaginary 
time formalism, which amounts to a compactification of the Euclidean time coordinate. The space-time becomes 
a topological cylinder. (As usual in this context, we consider only flat space-times without boundary.) Now, the 
heat kernel describes how an initial Dirac delta function in the space-time manifold spreads out as the proper time 
passes, with the Klein-Gordon operator K acting as a Laplacian operator. As it is known, the standard small r 
asymptotic expansion is insensitive to global properties of the space-time manifold. This means that the space- 
time compactification, and hence the temperature, will not be seen in the strict expansion in powers of r. (As a 
consequence, the ultraviolet sector and hence the renormalization properties of the theor y an d the quantum anomalies 
are temperature independent, a well-known fact in finite temperature field theory p3. |25|.1 Within a path integral 
formulation of the propagation in proper time, this corresponds to an exponential suppression (namely, of order 
e -/3 /4t^ c f ||2.5| l) of closed paths which wind around the space-time cylinder. The compactification is made manifest 
if instead of counting powers of r one classifies the contributions by their mass dimension. The corresponding thermal 
Seeley-DeWitt coefficients will then be powers of r but with exponentially suppressed r-dependent corrections. As 
a result of the compactification, the new expansion will not be Lorentz invariant, although rotational invariance will 
be maintained. In addition, we find coefficients of half-integer order which at zero temperature can appear only for 
manifolds with boundary (as distributions with support on the boundary 26]). Such half-order terms vanish in a 
strict proper time expansion. 

Another relevant issue is the preservation of gauge invariance. At zero temperature the only local gauge covariant 
quantities available are the matter fields, the field strength tensor and their covariant derivatives. However, at finite 
temperature there is a further gauge covariant quantity which plays a role, namely, the (untraced) thermal Wilson 
line or Polyakov loop. Since temperature effects in the imaginary time formalism come from the winding around the 
space-time cylinder, the Polyakov loop appears naturally in the thermal heat kernel. Our calculation, anticipated in 
27], shows that the thermal heat kernel coefficients at a point x become functions of the untraced Polyakov loop that 
starts and ends at x. Although such a dependence is consistent with gauge invariance at finite temperature, is not 
required by it either. Nevertheless, there is a simple argument which shows that the heat kernel expansion cannot be 
simply given by a sum of gauge covariant local operators (albeit with Lorentz symmetry broken down to rotational 
symmetry). For the Klein-Gordon operator describing a gas of identical particles free from any external fields other 
than a chemical potential (plus a possible mass term) it is obvious that such chemical potential (which can be regarded 
as a constant c-number scalar potential Aq) has no effect through the covariant derivatives, and so it is invisible in 
the gauge covariant local operators. However, it is visible in the Polyakov loop, and it is only in this way that the 
effective action, or the grand-canonical potential, and hence the particle density, can depend on the chemical potential. 
The dependence of the thermal heat kernel coefficients on the Polyakov loop was overlooked in previous calculations 
p8l l29l ]. although was made manifest in particular cases and configurations in [3(j. Of course, the relevance of the 
Polyakov loop is well-known in quarkless QCD at high temperature, where it is the order parameter signaling the 
presence of a deconfining phase j5| . The determination of the effective action of the Polyakov loop after integration 
of all others degrees of freedom has been pursued e.g. in |3l|]. Our results imply that, because the formulas are quite 
general and should hold for any gauge group, the Polyakov loop must be accounted for, not only in the color degrees of 
freedom and at high temperature, but also in other cases such as the chiral flavor group with vector and axial-vector 
couplings and at any finite temperature |l4j . The thermal heat kernel expansion is derived in Section ITT1 

In Section HTTI we apply the previous technique to the computation of the effective action of QCD at finite temperature 
to one-loop. Here we refer to the effective action in the technical sense of generating function of one-particle irreducible 
diagrams. For the quark sector (we consider massless quarks for simplicity) the method applies directly by taking 
as Klein-Gordon operator the square of the Dirac operator and using an integral representation for the fermionic 
determinant. In the gluon sector, the fluctuation operator is of the Klein-Gordon type in the Feynman gauge, 
and so the technique applies too, but this time in the adjoint representation of the gauge group and including the 
ghost determinant. The calculation is carried out using the covariant background field method. To treat ultraviolet 
divergences dimensional regularization is applied, plus the MS scheme. We have also made the calculation using 
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the Pauli-Villars scheme as a check. In this computation the background gauge fields are not stationary, and this 
allows to write expressions which are manifestly invariant under all gauge transformations (recall that in the time 
compactified space-time there are topologically large gauge transformations The result is expressed using gauge 

invariant local operators, including operators of up to dimension 6, and the Polyakov loop Q(x). This is done for 
arbitrary SU(A) (N being the number of colors). For SU(2) and SU(3) the traces on the color group are worked out, 
to dimension six for SU(2) and to dimension four for SU(3). In our expansion the dependence on the Polyakov loop 
is treated exactly (we keep all orders in an expansion in powers of log(f2)) but the expansion in covariant derivatives 
is truncated. In particular the time covariant derivative is not kept to all orders and in this we differ from .33j (who 
only treat SU(2) for stationary configurations). 

As it is known, the effective action of perturbative QCD at finite temperature contains infrared divergences due 
to the massless gluons in the chromomagnetic sector |34l l35j . Such divergences come from stationary quantum 
fluctuations which are light even at high temperature, whereas the non stationary modes become heavy, with an 
effective mass of the order of the temperature T, from the Matusbara frequency. So, the procedure which has been 
devised to avoid the infrared problem is to integrate out the heavy, non-stationary modes to yield the action of an 
effective theory for the stationary modes, i.e. of gluons in three Euclidean dimensions [Tol fTTL I3H I37I l38l I3H Eft l4ll |. 
In this way one obtains a dimensionally reduced theory C%d- (One can go further and integrate out the chromoelectric 
gluons which become massive through the Debye mechanism. We do not consider such further reduction here.) By 
construction, £30 reproduces the static Green functions of the four-dimensional theory C^d- Of course, the infrared 
divergences will reappear now if this action is used in perturbation theory. However, living in a lower dimension, C 3 jj 
is better behaved in the ultraviolet and also more amenable to non perturbative techniques, such as lattice gauge 
theory. The parameters of C^o (masses, coupling constants) can be computed in standard perturbative QCD since 
they are infrared finite, coming from integration of the heavy non stationary modes, although they are scale dependent 
due to the standard ultraviolet divergences of four-dimensional QCD. Section Hvl is devoted to obtaining the action of 
the reduced theory. This is easily done from the calculation of the effective action in Section lTTTl bv removing the static 
Matsubara mode in the gluonic loop integrations. This theory inherits the gauge invariance under stationary gauge 
transformations of the four-dimensional theory, but a larger gauge invariance is no longer an issue since more general 
gauge transformations would not preserve the stationarity of the fields. In addition, at high temperature fluctuations 
of the Polyakov loop far from unity (or from a center of the gauge group element in the quarkless case) are suppressed 
and so it is natural to expand the action in powers of Aq. We obtain the action up to operators of dimension six 
included (counting each gluon field as mass dimension one) and compare with existing calculations to the same order 
quoted in the literature [H El El El El • The relevant scales K T M E for the running coupling constant in the 
high temperature regime are identified and reproduced |43j |. For the dimension six terms, in the gluon sector we find 
agreement wit h | 42l if the Polyakov loop is expanded in perturbation theory and in the quark sector we reproduce 
the results of |44| for the particular case considered there (no chromomagnetic gluons and no more that two spatial 
derivatives). We give the general result for SU(iV) and simpler expressions for the cases of SU(2) and SU(3). 

The heat kernel and the QCD parts of the paper may interest different audiences, the first one being more method- 
ological and the second one more phenomenological, and to some extent they can be read independently. The QCD 
part does not require all the details of the derivation of the thermal heat kernel expansion but only the final formulas. 
In fact, one of the points of this paper is that the thermal coefficients need not be computed each time for each 
application but only once, and then applied in a variety of situations. 

II. THE HEAT KERNEL EXPANSION AT FINITE TEMPERATURE 

A. The Polyakov loop and the heat kernel 

We will consider Klein-Gordon operators of the form 

K = M(x)-Dl, D tl = d^ + A^(x). (2.1) 

M(x) is a scalar field which is a Hermitian matrix in internal space (gauge group space), the gauge fields A^(x) are 
antihermitian matrices. K acts on the particle wave function in d + 1 Euclidean dimensions and in the fundamental 
representation of the gauge group. At finite temperature in the imaginary time formalism the time coordinate is 
compactified to a circle, i.e., the space-time has topology Aid+i = S 1 x Md- Correspondingly, the wave functions 
are periodic in the bosonic case, with period f3 (the inverse temperature) antiperiodic in the fcrmionic case, and the 
external fields M, A^ are periodic. 

In order to obtain the heat kernel (x\e~ rK \x) (a matrix in internal space) we use the symbols method, extended to 
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finite temperature in |45ll46j: For an operator / = f{M, D M ) constructed out of M and D^, 

(x\f(M,D IM )\x) = ~'£ J ^<z|/(M,D M +ifv)|0>. (2.2) 

Here po are the Matsubara frequencies, 2im/(3 for bosons and 1-n{n + \)/ ft for fermions, and the sum extends to 
all integers n. On the other hand, |0) is the zero momentum wave function, (ir|0) = 1. The matrix valued function 
(x\f(M, + zp M )|0) is the symbol of /. Is important to note that this wave function is periodic (in fact constant) 
and not antiperiodic, even for fermions. The antiperiodicity of the fermionic wave function is only reflected in the 
Matsubara frequencies in this formalism. Whenever the symbols method is used, <9 M acts on the periodic external 
fields. Ultimately acts on |0) giving zero (this means in practice a right-acting derivative operator). 

In order to introduce the necessary concepts gradually, and to provide the rationale for the occurrence of the 
Polyakov loop in the simplest case, in what remains of this subsection we will consider the case of no vector potential, 
space independent scalar potential and constant c-number mass term, 

A(x) = 0, A Q = A (x ), M(x)=m 2 , [m 2 , ]=0. (2.3) 

This choice avoids complications coming from the spatial covariant derivatives and commutators at this point of the 
discussion. The result will be the zeroth order term of an expansion in the number of commutators [D^, ] and [M, ]. 
An application of the symbols method yields in this case 
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- (47rr w 2/3 ^>|e^° + ^ 10} • (2.4) 

(After the replacement D — > D + p dictated by 11'. 21) Di = di can be set to zero due to |0).) 

The sum over the Matsubara frequencies implies that the operator ^ ~52 Po e T ( Do+m) ' >2 is a periodic function of Do 

with period 2iri/f3, thus it is actually a one- valued function of e~P D ° . This can be made explicit by using Poisson's 
summation formula, which yields 



R ^ 



^ , o w = _L_ D±)W ^ (2 . 5) 

po K J feez 

(± for bosons or fermions, respectively). This observation allows to apply the operator identity |46j 

e f3do e -PD = q^j : (2 6) 

where ft(x) is the thermal Wilson line or untraced Polyakov loop 

/ rxo+P \ 

n(x) = Texp - / A (x' ,x)dx' (2.7) 



(T refers to temporal ordering and the definition is given for a general scalar potential Aq(x).) The Polyakov 
loop appears here as the phase difference between gauge covariant and non covariant time translations around the 
compactified Euclidean time. Physically, the Polyakov loop can be interpreted as the propagator of heavy particle 
in the gauge field background. The identity (|2.6|l is trivial if one chooses a gauge in which Aq is time independent 
(which always exists globally) since in such a gauge Q — e~P A °, and Do, Ao and do all commute. The identity itself 
is gauge covariant and holds in any gauge |46|. 

The point of using (|2.6|) , is that the translation operator in Euclidean time, e@ da , has no other effect than moving 
Xo to xo + and this operation is the identity in the compactified time, 

e 1390 = 1 , (2.8) 

(even in the fermionic case, recall that after applying the method of symbols the derivatives act on the external fields 
and not on the particle wave functions), so one obtains the remarkable result 



-PD 



fl(x). (2.9) 
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That is, whenever the differential operator Do appears periodically (with period 2m/ 0), it can be replaced by the 
multiplicative operator (i.e. the ordinary function) — ^ \og(Vl(x)). The many-valuation of the logarithm is not effective 

due to the assumed periodic dependence. Another point to note is that D (or any function of it), acts as a gauge 
covariant operator on the external fields F(xo,x), and so transforms according to the local gauge transformation at 
the point (xq,x). Correspondingly, the Polyakov loop, which is also gauge covariant, starts at time £0 and not at 
time zero in (|2.7|) : this difference would be irrelevant for the traced Polyakov loop, but not in the present context. 
An application of the rule (|2.9|l . yields in particular, 



1 

Po V ' k£? 



e r(D 0+lpo f = f_ y ( ± )fe n fc e -fe a /3 2 /4r _ (2 1Q) 



More generally, 

E + A>) = £ /(ipo - ^ log(fl)) , (2.11) 

Po Po 

provided the sum is absolutely convergent, so that the sum is a periodic function of Dq. Thus it will prove useful to 
introduce the quantity Q defined as 

Q = i Po + D = i Po - \ log(fi) . (2.12) 

P 

The second equality holds in expressions of the form (|2.11|l . (Note that the two definitions of Q are not equivalent in 
other contexts, e.g. in J2 po fi (Q)Xf 2 (Q) unless [£>o,X] = 0.) 
The heat kernel in (|2.4|) becomes 

x ' Po 

' - Tm Vo(«). (2.14) 



( 47rr )(d+l)/2 

In the first equality we have removed the bracket (x\ ■ |0) since for multiplicative operators like il(x), this brackets 
just picks up the value of the function at x. In the last equality we have used the definition of the functions (/? n (f2) 
which will appear frequently below 

Vn {n-r/(3 2 ) = (H^l^W^ 2 , Q^po^log(fi). (2.15) 

Po 

Note that there is a bosonic and a fermionic version of each such function, and the two versions are related by the 
replacement ft — > — fl. As indicated, these functions depend only of the combination t / 1 . In the zero temperature 
limit, the sum over po becomes a Gaussian integral, yielding 

^.») = {<-*>""<" -D" ( , 16) 

As can be seen for instance from (|2.1U|) . in this limit only the k = mode remains, whereas the other modes become 
exponentially suppressed, either at low temperature or low proper time r. 

The result in (|2.14() is sufficient to derive the grand canonical potential of a gas of relativistic free particles. For 
definiteness we consider the bosonic case 01 • The effective action (related to the grand canonical potential through 
W — flflgc) is obtained as 

r°° fj T 

W = Trlog(iT) = -Tr / — (x\e~ TK \x) (2.17) 
Jo T 

K includes a chemical potential Ao — —i/J, as unique external field, a the corresponding Polyakov loop is O = exp(z/3/z). 
Using (|2.14() . subtracting the zero temperature part (which corresponds to setting ipo — > 1) and carrying out the 
integrations yields the standard result [24J 

W = N J ^0^. [bg (l - e -f^-^ + bg (l - e -?K+"))] . (2.18) 
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N is the number of species and lu^ = \/k 2 + m 2 . 

In next subsection, after the introduction of more general external fields, we will consider expansions in the number 
of spatial covariant derivatives and mass terms. At zero temperature, the derivative expansion involves temporal 
derivatives as well, as demanded by Lorentz invariance, but such an expansion is more subtle at finite temperature. 
The direct method would be to expand in powers of Dq in l|2.4|l . however, this procedure spoils gauge invariance (e.g. 
Dq|0) = Ao|0) is not gauge covariant). As a rule, giving up the periodic dependence in Do breaks gauge invariance 
|46j . One can try to first fix the gauge so that Aq is stationary, and then expand in powers of Ao . This is equivalent 
to expanding in powers of log(fi). By construction this procedure preserves invariance under infinitesimal (or more 
generally, topologically small) gauge transformations, however, it does not preserve invariance under discrete gauge 
transformations (0,E3 and 1111 l51 below). This is because log(fi) is many- valued under such transformations. An 
expansion in the number of temporal covariant derivatives which does not spoil one-valuation nor gauge invariance is 
described next. 



B. The diagonal thermal heat kernel coefficients 

Here we will consider the heat kernel expansion at finite temperature in the completely general case of non trivial 
and non Abelian gauge and mass term fields, A fl (x), M(x). 

First of all one has to specify the counting of the expansion. At zero temperature, the expansion is defined as 
one of (x\e~ rK \x) in powers of r [after extracting the geometrical factor (47rT)~( d+1 )/ 2 ]. Each power of r is tied to 
a local operator constructed with the covariant derivatives D M and M{x) (cf. I|2.23[l and l |2.24| l). The heat kernel 
e~ rK is dimensionless by assigning engineering mass dimensions —2, +1, +2 to r, and M, respectively. So, at zero 
temperature, the expansion in powers of r is equivalent to counting the mass dimension carried by the local operators. 

At finite temperature there is a further dimensional quantity, /3, the two countings are not longer equivalent 
and one has to specify the concrete expansion to be used. It is well-known that the finite temperature corrections 
are negligible in the ultraviolet region, so that for instance, the temperature does not modify the renormalization 
properties of a quantum field theory [3 E3 and also the quantum anomalies are not affected 0, • The ultraviolet 
limit corresponds to the small r limit in the heat kernel. As noted before, and can be seen e.g. in i|2.10[l . the finite j3 
and small r corrections are of the order of / 4r or less, and so they are exponentially suppressed. Of course, the 
same exponential suppression applies to the low temperature and finite t limit. This implies that a strict expansion 
of the heat kernel in powers of r will yield precisely the same asymptotic expansion as at zero temperature. In order 
to pick up non trivial finite temperature corrections we arrange our expansion according to the mass dimension of 
the local operators. In this counting we take the Polyakov loop fl, and M as zeroth, first and second order, 
respectively. In addition one has to specify that fl(x) is at the left in all terms (equivalently, one could define a similar 
expansion with Q always at the right). This is required because the commutator of f2 with other quantities generates 
commutators [Do, ] which are dimensionful in our counting. After these specifications the expansion of (x\e~ rK \x) for 
a generic gauge group is unique and well-defined and full gauge invariance is manifest at each order. In this expansion 
the terms are ordered by powers of r but with coefficients which depend on 1 jr and fl, 

{x\e-^ M - D ^ \x) = (47rr)"( £i + 1 )/ 2 £ a T n {x)r n . (2.19) 

n 

From the definition it is clear that the zeroth order term for a general configuration is just 

a£(x) = <po(n(x);T/f3 2 ) (2.20) 

already computed in the previous subsection (cf. I|2.14|I % 1. This is because when the particular case (|2.3() is inserted in 
the full expansion all terms of higher order, with one or more [D^, ] or m 2 , vanish identically. 

For subsequent reference we introduce the following notation. The field strength tensor is defined as F M „ = [D^, D u \ 
and likewise, the electric field is Ei = Foi- In addition, the notation means the operation [D^, ]. Finally we will 
use a notation of the type X^ ua to mean D ll D v D a X = [D^, [D v , [D a , X]]], e.g. M 00 = D%M, F nl ,„ = D a F^ v . 

The method for expanding a generic function (x\f(M, D^x) has been explained in detail in [ifj. We have applied 
this procedure to compute the heat kernel coefficients to mass dimension 6. However, for the heat kernel there is an 
alternative approach which uses the well-known Seeley-DeWitt coefficients at zero temperature. This is the method 
that we explain in detail here. The idea is as follows. The symbols method formula 12.21) is applied to the temporal 
dimension only 

(x\e-« M - D V\x) = ~Y l (x\B^ M - t > a - D n\x), Q = i Po +D (2.21) 

Po 
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(The bracket (xq\ |0), associated to the Hilbert space over xq, is understood although not written explicitly.) This 
implies that we can use the standard zero temperature expansion for the d-dimensional heat kernel with effective 
Klein-Gordon operator 

K = Y- D 2 , Y = M -Q 2 . (2.22) 

In this context Y is the non Abelian mass term, because, although it contains temporal derivatives (in Q), it does 
not contain spatial derivatives and so acts multiplicatively on the spatial Hilbert space. The standard heat kernel 
expansion gives then 

oo 

(x\e-^ Y - D ^\x) = (4nT)- d ^ J2 a n(Y, A)r" , (2.23) 

71=0 

where the coefficients a n (Y,Di) are polynomials of dimension 2n made out of Y and Di — [Di, ]. To lowest orders 
dill 

00 = 1 , 

01 = -Y, 

a 2 = -Y 2 - -Y it + —F 2 , 
2 6 12 13 

as = -g^ 3 + j^{Y, Ya} + —Y 2 - —Y Hjj - — [Fuj , Yj] 

-™{ Y > F «) ~ Kn Fi 3 YF v + JS F ijk ~ ^ F v F jkF kt + Yon F iij + fin^'' Fkki i } ■ ( 2 - 24 ) 



(As noted before Yu = D 2 Y, F tjk = DiFj k , etc.) 

Eq. (|2.23|l inserted in (|2.21|) is of course correct but not very useful as it stands. For instance, for the zeroth 
order, the expansion in (|2.23|l would be needed to all orders to reproduce the simple result H2.2()(l . since e T ® is not a 
polynomial in Q. In view of this, we consider instead 

oo 

(x\e-^ M ~Q 2 - D ^ \x) = (^r)-^ 2 e TQ2 ~a n (Q 2 ,M, A)r" , (2.25) 

n=0 

which introduces a new set of polynomial coefficients a n (Q 2 , M, Di). By their definition, it is clear that these coef- 
ficients are unchanged if "Q 2 " is everywhere replaced by "Q 2 + c-number" . This implies that in a n the quantity Q 2 
appears only in the form [Q 2 , ]. This is an essential improvement over the original coefficients a n , since each [Q 2 , ] 
will yield at least one _D , and so higher orders in [Q 2 , ] appear only at higher orders in the heat kernel expansion 1 . 
The calculation of the coefficients a n (Q 2 , M, Di) follows easily from the relation 

oo oo 

]T a n r n = e^ 2 ]T ~a n r n . (2.26) 

n=0 n=0 

If one takes the expression on the l.h.s. and moves all Q 2 blocks to the left using the commutator [Q 2 , ], two types 
of terms will be generated, (i) terms with Q 2 only inside commutators and (ii) terms with one or more Q 2 blocks at 
the left. The terms of type (i) are those corresponding to J2 n ^ T "- To lowest orders one finds 

d = 1 , 
Si = —M , 

a 2 = hi 2 - hl u + ±F 2 + 1 [Q 2 ,M] + ^(Q 2 ) lt . (2.27) 

Once the d n coefficients are so constructed one has to proceed to rearrange l|2.25|l as an expansion in powers of M, 
Di, and Dq. The expansions in M and Di are already inherited from l|2.23() . It remains to expand [Q 2 , ] in terms of 



1 This kind of resummations are standard also at zero temperature to move e.g. the mass term e tM to the left and leave only a [M, ] 
dependence in the coefficients Il7l . 
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[Q, ], or equivalently, in terms of Dq — [Do, ] since the quantities Q and Dq differ by a c-number. To do this, in the 
a n coefficients Q is to be moved to left, introducing Dq, until all the terms so generated are local operators made out 
of Du and M, and all uncommutated Q's are at the left. E.g. 

a 2 = ^M 2 -^Mu + ^F 2 -^M o + ^E 2 + ^E OH + QM -^QE u . (2.28) 

(Recall that Ei stands for the electric field Foi.) We can see two types of contributions in 62, namely, those without 
a Q at the left and those with one. If Q is assigned an engineering dimension of mass, all the terms are of the same 
dimension, mass to the fourth. However, in our counting only the dimension carried by and M is computed, and so 
the two types of terms are of different order, namely, mass to the fourth and mass to the third, respectively. Indeed, 
when 02 is introduced in (|2.25[) (i.e. it gets multiplied by e T< ^') and then in (|2.21|) (the sum over the Matsubara 
frequencies is carried out) we will obtain the contributions (using X^ Po Q n e T ® ~ ip n ) 

«2 - po(fl) Qm 2 - ~M U + - ^Moo + ^E 2 + ^E 0l ^j r 2 + <^(Q) \M - ^E tl ^j r 3 / 2 . (2.29) 

These are contributions to the thermal heat kernel coefficients and a,y 2 , respectively, introduced in i|2.19[l . Note 
the presence of half-integer order coefficients from terms with an odd number of Q's. 

As we have just shown, each zero temperature heat kernel coefficient in (|2.23|1 allows to obtain a corresponding 
coefficient a,k with the same engineering dimension 2k. Such coefficient in turn contributes, in general, to several heat 
thermal coefficients (with mass dimension 2n). Let us discuss in detail to which contributes each The 
change from engineering to real dimension comes about because some terms in dk contain factors of Q at the left 
which do no act as -Do and so count as dimensionlcss. Therefore it is clear that for given k, the allowed n satisfy 
n < k, the equal sign corresponding to terms having all Q's m commutators. On the other hand, the maximum 
number of [Q 2 , ]'s in dk (k > 0) is k — 1, and from these, at most k — 1 uncommutated Q's can reach the left of the 
term. This yields the further condition k < 2n — 1. Note further that a factor Q e gives rise to a coefficient <pe(£l) 
in a^. In summary, in the computation of the thermal coefficients a„ up to n = 3 (mass dimension 6), we find the 
following scheme 



a r 


- CL r 


- ip a% 








Ol - 




- cp of 








a 2 - 






h (pi a$ /% 






a 3 - 


- a 3 - 


- ^4 - 


I" <Pl a 5/2 ~ 








- 04 r 


- Lp aj - 








a 5 - 


" a 5 - 




h </?i a^ /2 - 


- 04 - 


h <^3 a 7 / 2 



<^ 4 (2.30) 

The mixing of terms is a nuisance that does not occur at zero temperature, however, it cannot be avoided: Q contains 
Po and must count as zeroth order (otherwise, if Q were of order one the expansion would consist of polynomials in Q 
and the sum over po would not converge) . On the other hand, counting po as zeroth order and Do as first order even 
when it is inside Q results in a breaking of gauge invariance, as we noted at the end of the previous subsection. The 
fact that Q counts as dimensionlcss and Do as dimension 1 is necessary to have an order by order gauge invariant 
expansion. This counting is well defined provided that all fi's are at the left (for instance) of the local operators (cf. 
(|2.36[) and discussion below). 

From l|2.30|l we can see that we do not need the complete zero temperature coefficients 04 and 05. requires only 
terms Y n , with n = 2, 3, 4 in 04 (Y, Di) and n = 4, 5 in as(Y, Di). We have extracted the zero temperature coefficients 
from [l8|. These authors actually provide the traced coefficients b n (x) defined by 

00 

Tr (e~ T ^ Y - D ^ = (4irT)- d / 2 J2 / d d xtv{b n )r n , (2.31) 

n=0 J 

where Tr is the trace in the full Hilbert space of wave functions and tr is the trace over the internal space only. The 
coefficient a n is obtained by means of a first order variation of b n +i (cf. (|2.41l) '). The advantage of this procedure is 
that the traced coefficients are much more compact and better checked. 

As we have said, we have computed the thermal heat kernel coefficients up to and including mass dimension 6 by 
the procedure just described and also by that detailed in [4^. This latter approach uses the symbols method for 
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space and time coordinates and so computes the coefficients from scratch (in passing it yields the zero temperature 
coefficients as well). We have verified that the two computations give identical results after using the appropriate 
Bianchi identities (in practice the method of tends to give somewhat more compact expressions). The results are 



as follows 




T 
u 




T 
a l/2 


= o, 


T 
«1 


= -<p M, 


T 

a 3/2 


= <pi \ Mo - 


T 
«2 


T=0 

= tp a 2 - 


T 
a 5/2 


= | (2^i + 



—Ej. 
3 

\y 2 {E% + E 0ii - 2M 00 ) , 
6 



l^Mou - \<p x (2M M + MM ) (2.32) 
6 3 



~<Pi ({Mi, Ei} + {A/, E u }) - \ hp x + -v? 3 J E 00ii - 





2 \ 






(^1- 


- -tp 3 J E 0i Ei - 




-IV 3 ) 



?0i + 2Q<Pl[Ej>Fiij] 



tp aJ-° - ( -Tp 2 - -^<p 4 J Moooo ~ (3M oii - 15M 00 M - 5MM 00 - 1 5A/ 2 

+4{Af, E|} + 2^ 4 M^ + AMEoa + &E Qii M + 4M i £ CH + 6E Ql M l 
+7M E U + 3E U M + 6M 0i Ei + 4E t M 0l 



3^ 2 - iV 4 ) ^ 00< + (sl^ 2 - ^ ' ^ 



IgO^ 2 ^^*' + 4{_Ei, Sijj-} + 5i? 2 + 4_E 2 + AFoujEj — 2EjFouj — 2E ni jFij 

[Eij, Foij] — AEoiFjji + 2FjjiEoi + 2EiFijEj + 2{EiEj, Fij} + 7FooijFij + 3FijFnoij + SF^j 



In these formulas a^ =0 stands for the zero temperature coefficient. These are the same as those in (|2.24() but using 
M instead of Y and space-time indices instead of space indices, e.g. a 2 =0 — \M 2 — \M^ + j^F^. For convenience 
we have introduced the auxiliary functions 

4 (— 2)™ 

Tp 2 = if Q + 2ip 2 , Tp 4 = ip - -ip 4 , Tp 2n = (f - _ jTjy l2n , (2.33) 

which vanish at r//3 2 =0. Due to the Bianchi identity there is some ambiguity in writing the terms. We have chosen 
to order the derivatives so that all spatial derivatives are done first and the temporal derivatives are the outer ones. 
This choice appears naturally in our approach and in addition is optimal to obtain the traced coefficients 6^ since the 
zeroth derivative of the Polyakov loop vanishes (cf. (|2.3t)|) below), and so terms of the form ip n Xo do not contribute 
to the traced coefficients upon using integration by parts. The terms Oq, af, a^, 2 and a 2 were given in p7j . 



C. The traced thermal heat kernel coefficients 



The zero temperature traced heat kernel coefficients have been introduced in (|2.3f() (for the ci-dimensional operator 
Y — Df). Of course, the choice b n = a n would suffice, however, exploiting the trace cyclic property and integration 
by parts more compact choices are possible. At lowest orders the coefficients can be taken as (we give the formulas 
for K = M — D 2 at zero temperature; the heat kernel coefficients are dimension independent) \li 

bo = f, 
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h = —M , 

b 2 = -M 2 + —F?„. , (2.34) 
1 2 12 * v 



b% — — M- — — — — F tn ,AlF in , — — F,,,,,, ~t~ — Fm.FjjnFrv 
■ q 12 12 60 90 



By construction a n — b n is a commutator which vanishes inside Tr. Likewise, we can introduce the traced coefficients 
at finite temperature, 



Tr 



with b„ simpler than a^. Once again we choose a canonical form for these coefficients where a function of f2 put at the 
left is multiplied by a local operator (i.e. an operator made out of M and D^). To simplify the traced coefficients and 
bring them to the canonical form we need to work out the commutators of the form [X, f(Q)] (in particular D^f(fl)) 
as a combination of terms of the type function of f2 times local operator. As shown in Appendix [S] the rules are 
as follows: let / denote a function of £1 (e.g. </?„(0)) and let be its n-th derivative with respect to the variable 
-log(0)//3, then 

D f = 0, 

D t f = -f'Ei + \f"E 0i - -fWEaoi + • • • , (2.36) 



[XJ] = -f'X + lf"X 0Q -lf^X 



2 J J1 uu 2 1 J J1 -000 

These formulas imply that, unlike the zero temperature case, the cyclic property mixes terms of different order at 
finite temperature. This is because, as noted above, D has dimensions of mass whereas fl counts as dimensionless. 
So for instance, (fo(£l) is of order zero and Di is of first order, yet Diipo(Q) contains terms of all orders, starting with 
dimension 2. As we will discuss below, this implies that there is a certain amount of freedom in the choice of the 
traced coefficients. To apply these commutation rules to we further need the relation 

f'n = Vri^n+i + rup n -i) . (2.37) 

Using these rules we can apply integration by parts and cyclic property to the previously computed coefficients 
and choose a more compact form for them valid inside the trace. In this way we obtain, up to mass dimension 6, 

bo = fo, 

bl/2 = 0, 

b\ = -<poM, 

bI/2 = 0, (2.38) 

b 2 = f b 2 - -^P2 E i > 

65/a = -\<Pi{Mi,Ei} , 

bj = foh + \<fr U^ 2 + E^IE, + ^El + ^F^ - ^ /. , f , .. /.\ ) + (-^ 4 - j^ 2 J E 2 0l . 

This is the main result of these Section, where the <p n functions are given in l|2.15|l and (|2.33|) . In these formulas the 
b n are the zero temperature coefficients given in (|2.34|l . We note that the coefficient bj above is not identical to that 
given in ^fy. (The coefficient in [23 corresponds to replace (^0^3 above by ifiob' 3 , where 63 differs from 63 in (|2.34|l by 
a cyclic permutation.) The two versions of bj differ by higher order terms. In what follows we use the coefficient in 
fTM . 

Several remarks should be made about these expressions. Either at zero or finite temperature there is an ambiguity 
in the choice of the traced coefficients b^, however, the ambiguity is essentially larger at finite temperature. Indeed, 
writing the expansion as 

Tr U-< M -K^ = (47rr)-( d+1 )/ 2 £ B T n r n , (2.39) 
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we find that, although b n is ambiguous, B^= a is not. This is because at zero temperature the expansion is tied to a 
series expansion in powers of a parameter (say, r). At finite temperature the expansion is not tied to a parameter 
(it is rather a commutator expansion) and so the ambiguity exits not only for 6^ but also for B^- For instance, bj 
above has been expressed in terms of the coefficient bi given in l|2.34[l . Nothing changes at zero temperature if we add 
to 62 since the addition is a pure commutator, however, in b 2 it would mean to add (poM^ which is no longer 
a pure commutator, thereby changing the functional B 2 . In fact, tpoM^^, which is formally of dimension 4, can be 
expressed as a sum of terms of dimension 5 and higher, using integration by parts and the commutation rules l|2.36[) . 
So the concrete choice of b 2 affects the form of the higher orders, b^ 2 , bj, etc. 

Taking into account this ambiguity, our criterion for choosing the traced coefficients has been to recursively bring 
the &Jj to a compact form. We observe that inside the trace (upon the applying the commutation rules) a,y 2 is a sum 

'^j 2 — *~>. J- Uj 2 i a.u.gxj.i.cuijc<a vvilii inu t^i-Liio £,c±icicLLcd from ^3^2' 

T 

5/2' 



of terms of dimension 4 and higher, so we choose b^ 2 = 0. Then a 2 , augmented with the terms generated from a T 
is brought to the most compact form. This in turn produces higher order terms which are added to aT, 2 , and so on 



Of course, this is not the only possibility, since taken a b^ to be simplest may imply a greater complication in the 
higher order coefficients. For instance, as can be shown, it is possible to arrange the expansion so that all half-order 
traced coefficients vanish. E.g. b^ 2 can be removed at the cost of complicating bj . 

It should be clear that the ambiguity in the expansion in l|2.39[l does not affect its sum but only amounts to a 
reorganization of the series. On the other hand the untraced coefficients are not ambiguous: once brought to their 
canonical form they are unique functionals of M and A^. 

The heat kernel is symmetric under transposition of operators, the b^ have been chosen so that this mirror symmetry 
holds at each order. 

As is well-known ^3] not only the allow to obtain the 6j but also the converse is true. By their definition 

{ X \e-^- D D\x) = -'-j^Tr ( e -*<"-*S)) . (2.40) 
Using the expansions in both sides, one finds at zero temperature (using l|2.39|l 'l 



a 



r =o (a . ) = _^4. (2.41) 
SM(x) y ' 



5B. 



At finite temperature, the variation of b\ contributes not only to a r [_ 1 but also to all higher order coefficients, in 
general. So we have instead 

«*<•> = -«k> E BlrM - <2 ' 42) 

where in the r.h.s. only the terms of dimension 2n are to be retained and k takes integer as well as half-integer values. 
We have checked our results by verifying that this relation holds for our coefficients. 



III. THE ONE-LOOP EFFECTIVE ACTION OF CHIRAL QCD AT HIGH TEMPERATURE 

Here we will apply the thermal heat kernel expansion just derived to obtain the one-loop effective action of QCD 
with massless quarks in the high temperature region. We remark that the effective action we are referring to is the 
standard one in quantum field theory, namely, the classical generator of the one-particle irreducible diagrams. As a 
consequence our classical fields may be time-dependent. The quantum effective action in the sense of dimensional 
reduction [33, as an effective field theory for the static modes is of great relevance in high temperature QCD and is 
also discussed below, in Section ITY1 We will use the background field method, which preserves gauge invariance [50| . 
The Euclidean action is 

S = ~W I d4xtT{F ^ ) + J ^vPv- 

Dfi = +A fJ- , with and F^ v = [D^jD^] antihermitian matrices of dimension N. They belong to the fundamental 
representation of the Lie algebra of the gauge group SU(A). 2 



2 Our point of view will be that A M itself is the quantum field, independently of any particular choice of basis in su(N). So the coupling 
constant g is also independent of that choice. Due to gauge invariance is not renormalized. 
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A. Quark sector 



In this subsection we work out the quark contribution which is somewhat simpler than the gluon contribution. 
(The latter requires the use of the adjoint representation, introduction of ghost fields and treatment of the infrared 
divergences.) Upon functional integration of the quark fields, the partition function of the system picks up the 
following factor from the quark sector 

Z g [A] = Det(p) N f = Det(^ 2 )^/ 2 , (3.2) 

where Nf denotes the number of quark flavors. (As usual, we have squared the Dirac operator to obtain a Klein- 
Gordon operator.) The corresponding contribution to the effective action is (we use the convention Z — e~ r ^) 

T q [A] = -^Tr log(^ 2 ) = -L J ^-Trexp (r p 2 ) =: J <fix£ q (x) , (3.3) 



„ , N N f f°° dr n 2e „ nT . 



(3.4) 



In this formula the Dirac trace is included in the q and tr refers to color trace (in the fundamental representation). 
The ultraviolet divergences at t = are regulated using dimensional regularization, with the convention D = 4 — 2e. 
As is standard in dimensional regularization, the factor /x 2e is introduced in order to deal with an effective Lagrangian 
of mass dimension 4 rather than 4 — 2e. 

To apply our thermal heat kernel expansion we need only to identify the corresponding Klein-Gordon operator. We 
use 

7/^=7^! 7m lv = + , tr Dirac (l) = 4 . (3.5) 

The expression 

- f= -Dl - \a^F^ (3.6) 

identifies — ^(J^F^ as the (square) mass term M of the Klein-Gordon operator in this case. A direct application of 
(|2.38(l . shows that b\ and b^ 2 cannot contribute (they have a single M and this cancels due to the trace over Dirac 
space). The other coefficients give, to mass dimension 6 included, 

b Q, q = 4( A>, 

bl, q = -\{^Fl v +Tp 2 E 2 ) , (3.7) 



uT _ _(,~- I?2 i — j?2\ 

3 

b L = ^ (%F^F„ X F X , + ^ - ±F*\ (±E 2 - ±F Q % - ^E^eA + ~ ] ^ 



In these formulas the functions cp n (defined in (|2.15(l and (|^23) correspond to their fermionic versions. All fields are 
in the fundamental representation. 

The required integrals over r in (|3.4|l are of the form 

f°° dr 

lU")-= -(4Vr) e rV r tM, H = l (3.8) 
Jo T 

where <p^ refers to the bosonic or fermionic version, respectively. In the quark sector the argument to will be the 
Polyakov loop in the fundamental representation, or, in practice, any of its eigenvalues. These integrals can be done 
in closed form (see Appendix [BJ). In particular 

W e2 ™) = (-1) W (^) 2£ {£) 21 T Jl±^±l±M [ C (l + 21 + 26, \ + v) + C(l + 21 + 2 C) \ v) 

-\<»<\- (3-9) 
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The integrals if n (u>) are one- valued functions of w, i.e. periodic in terms of v, however, to apply the explicit formula 

(|3.9|) v has to be taken in the interval — \ < v < |. The generalized Riemann ^-function C( z :<?) — J2^Lo( n + ( l)~ Z 
has only a single pole at z = 1 |5l|. so the dimensionally regulated integrals yield the standard pole of the type 1/e 
solely for the integrals lQ 2n , which appear in b^ q . 

We can now proceed to compute the contributions to the effective Lagrangian. The zeroth order requires IZ 2 o- 
Using the relation £(1 — n,q) = —B n (q)/n, n = 1, 2, . . ., with B n (q) the Bernoulli polynomial of order n [5jJ, one finds 

I-2.0 = - | ( j ) 4 + i/) + 0(e) (3.10) 

so the effective potential is 

C , q (x)=^N f T^(^-^tr[(l-4V 2 ) 2 ]^ , Q(x)=a 2 ™, ~\<V <\. (3.11) 

Here A^ is the number of colors, tr is taken in the fundamental representation of the gauge group and v is the matrix 
log(£!)/(27ri) with eigenvalues in the branch |77| < 1/2. This is the well known result 
The terms of mass dimension 4 have a pole at e = 0. Using the relation 

C(l+*,9) = --#?)+ 0(*)> (3-12) 
(where ip(q) is the digamma function) one finds 

Io,o = \ + log(47r) - 1E + 2 log(/i/3/47r) - +u)- -v) + 0(e) , 

T o,2 : = J o,o + 2I o,2 = -2 + 0(e) . (3.13) 

(For convenience, we have introduced the integrals if^- analogous to lf 2n m jUtE)) but using <p 2n instead of </?2n-) 

The terms (e -1 + log(47r) — je) in Iq come with ti(F^ v ) and are removed by adopting the MS scheme. We will 
discuss this in conjunction with gluon sector. After renormalization 

= ~lj^ N f tr [(21og(/i/47rT) - i>(\+V) - V>(| - V)) F 2 „ - 2E 2 ] . (3.14) 

Finally, the terms of mass dimension 6 in four space-time dimensions require I± , and I7-t- Using the relation 
i]j {n) (q) = (-l) n+1 n\Q(n+l,q) 51], one obtains 

^"(Jr) Wih + ^ + ^ih-^+Oit)' I- 2 = -21^,0 + 0(e), /- 3 = -4I 1 - + O(e). (3.15) 
(All these integrals are related through simple proportionality factors, as follows from i|B6|) ). This yields 

(V'(i+77)+<(!-77)) 



L3 - q{ } ~ (4^T2 tr 



1 , 1 , 1,1 

■712 I T?2 J^2 



X ' 45 ^^ v ^ vX ^ x ^ + 24^ MI/ 60 + 20^ OAtl/ 30^" + 15^^^ 



(3.16) 



In all these formulas v is the matrix log(f2)/(27ri) in the fundamental representation and in the branch \v\ < | in 
the eigenvalue sense. Note the hierarchy in powers of temperature, Cq ~ T 4 , £ 2 ~ T , £ 3 ~ T~ 2 , implying that the 
heat kernel expansion at finite temperature is essentially an expansion on k 2 /T 2 with k the typical gluon momentum. 
Terms of order T 2 are forbidden since there is no available gauge invariant operator of dimension two. 



B. Gluon sector 



In the background field approach [5(j the gluon field is split into a classical field plus a quantum fluctuation, i.e. 
— > + in the action (|3.1|l . As is standard in the effective action formalism, the appropriate currents are 
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added so that the classical field is a solution of the equations of motion (and so no terms linear in the fluctuation 
remain). The one-loop effective action corresponds then to neglect contributions beyond the quadratic terms in the 
quantum fluctuations and integrate over a^. (The quark fields are taken as pure fluctuation, so a M does not change 
the quark sector at one-loop.) 

The quadratic piece of the gluon action is 

S {2) = — a / d 4 xtr \-a v Dla v - 2a IM [F IMV ,a v ] - (D^) 2 ] . (3.17) 

Here all covariant derivatives are those associated to the classical gluon field . Note that the first two terms are of 
the standard Klein-Gordon form, but the last one is not. Before doing the functional integration over a M one has to 
fix the gauge of these fields. This implies adding a gauge fixing term and the corresponding Faddeev-Popov term jS^ 
in the action. We take the covariant Feynman gauge D^a^ = f(x), since the associated gauge fixing action precisely 
cancels the offending term (D^a^) 2 in l|3.17|l . After adding the ghost term one has 



S^ = -\ I d 4 xtr 
g 2 



-a u D 2 a v -2a^[F^,a v ]-CD 2 C . (3.18) 



The coupling constant has no effect here since it can be absorbed in the normalization of the fields. The ghost fields C 
and C are anticommuting (although periodic in Euclidean time) and are matrices in the fundamental representation 
ofsu(iV). 

The full effective action (to one loop) is 

= I d D xtv{Fl v )+T q [A]+Y g [A] , (3.19) 

where the first piece is the tree level action (accounting for renormalization; g is dimensionless), the second one is the 
quark contribution, obtained in the previous subsection, and the last term follows from functional integration over 
and C, C in HTWi . 

T g [A) = ^Trlog (-Dl ~ 2% v ) - Trlog =: J d A xC g {x) . (3.20) 

where = [£)„, ] and — [F^, }. From (|3.18|) . we can see that the Klein-Gordon operator over the gluon field 
a M acts on an internal space of dimension D x (N 2 — 1), where D = 4 — 2e is the number of gluon polarizations 
(including the two unphysical ones) and corresponds to the Lorentz index /i, and N 2 — 1 is the dimension of the 
adjoint representation of the group. _D„ and F^ v act in the adjoint representation. The covariant derivative of the 
Klein-Gordon operator is the identity in the Lorentz space whereas the "mass term" is a matrix in that space, namely, 
(M) M „ = —2F^ V . Similarly, the space of the Klein-Gordon operator over the ghost fields has dimension TV 2 — 1, the 
mass term is zero and the corresponding covariant derivative is just but in the adjoint representation. 
Applying once again the heat kernel representation, we have 

i r°° a t ,,2£ 

4,0°) = -o / "TI^E^O' (3.21) 



2j t (4ttt) d / 2 

where for convenience, the Lorentz trace over gluons as well as the ghost contribution are included in the coefficient 
g - tr denotes the color trace in the adjoint representation. A straightforward calculation yields 

bl g = (D-2) 9o (ti), 

bl 9 = (-2 + ^) - ^^(n)Ef , (3-22) 



,,4 D~2\~ 9 - l- 2 D-2~ 2 
°3, 3 - <M") ( ( 3 + qq J F ^ F v\Fxn + ^F XfiV ^— F^ v 



-Tp 2 {VL) ^-2F^ U + [El + F % - 2E i F ij E j 
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The coefficients bj g and bj, 2 g vanish, as do all terms with a single M, due to the Lorentz trace. The contributions with 
D — 2 come from pieces without M in (|2.38() and (|2.34|) . The effect of the ghost is to remove two gluon polarizations, 
D — > D — 2. Unlike the fermionic case, the thermal heat kernel coefficients depend explicitly on the space-time 
dimension through these polarization factors. In these formulas the functions tp n correspond to their bosonic versions. 
In addition, its argument f2 and all field strength tensor and covariant derivatives are in the adjoint representation. 

We can now proceed to the calculation of the effective Lagrangian. We note that the integrals over r are no 
different to those for the quark sector (see (|3.8|1 and Appendix 0, after the replacement v — > v — \ (coming from 
ifni^) = fn ( — Ll 0) an d so < v < 1 now: 



1 l.2n\ e I 



H0 



2tt 







2tt 



2f 



T(e + n + e+i) 



< v < 1. 



C(l + 21 + 2e, v) + ((1 + 21 + 2e, 1 - v) 



(3.23) 



In this way, for the effective potential one obtains 

7T 2 

£o, g (x) = yT 4 tr[5 4 (P) + B 4 (l-P)] ( 3 - 24 ) 

= -^T A (N 2 -l) + ^T 4 tr\v 2 (l-v) 2 ] , i? = log(fi)/(27rf) , 0<P<1. (3.25) 
45 3 L J 

This is also in agreement with the well known result 0. We emphasize that f2 and v are now in the adjoint 
representation as indicated by the notation tr. 

The mass dimension 4 piece of the effective Lagrangian, coming from 6^, requires I^q which is ultraviolet divergent 

and Jj"- which is UV finite (cf. (j3.13|l ). The finite pieces, in the MS scheme, are found to be 



(4tt) : 



■ti- 



ll 



g I 2 log( A1 /47rT) + ^j- - *l>(v) - ip(l -v)\ F* v - -E? 



< v < 1 



(3.26) 



On the other hand, the divergent contribution in the gluon sector, combined with that in the quark sector and the 
tree level Lagrangian yields (all terms have been multiplied by the factor /i 2<s to restore dimensions) 



1 



1 



2^ tr(i ^ ) ' N,ti- 



1 



log(47r) - 7b 



%tv(Fp - ^tr(F^) 



(3.27) 



Use of the SU(N) identity (|D5|I yields the renormalized tree level Lagrangian 

1 



C tIee (x)+£f v (x) + £f v (x) 



2g2( M ) 



tr(F 2 ) , 



(3.28) 



with the standard one-loop renormalization group improved in the MS scheme 



1 



1 

So 



1 



00 - + log(47r) - lE 



00 = 



1 



(4tt) 2 



11 2 
— TV — -N f 

q q J 



(3.29) 



guarantying the scale independence of (|3.27l) . Note that, due to gauge invariance, the classical fields do not 
need ultraviolet renormalization. (In the context of the dimensionally reduced effective theory, finite, temperature 
dependent, renormalization has been found to be useful in practice [33, |42|- See Section llVl ) 
Putting together all terms of mass dimension 4 (renormalized tree level plus one- loop), we find 



C 2 {x) 



1 



2g 2 M 
11 1 

12(4^)2 

1 1 



+ O log( A1 /4 7 rT) + 



1 1 



tr 



6 (4tt) 2 
(iP(D)+^(1-v))F 2 u 



3 (4tt) 2 
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N f tv[{i>{\+V)+^(\-V))FZ v } 



1 1 

-2 <V< 2 



<z?< 1. 



(3.30) 
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The terms of mass dimension 6 are easily obtained from the coefficient bj g and the integrals I^ , IT^ and 1^-, 

(^'(P)+^"(l-P)) 



11 1 ~ 



x ( J^hvF^F^ + -Fl^ v - —Fl^ + \Fo^ - + —EiF^Ej 



(3.31) 



45 3 30 w k u/xw 15 " 15 

Note again the hierarchy in powers of temperature, £o ~ T 4 , £2 ~ £3 ~ T~ 2 . 

C. Infrared divergence and other renormalization schemes 

The integrals may contain not only ultraviolet divergences but also infrared ones (corresponding to the large r 
region). Specifically, this happens if t > 0, n = and e 27rll/ = ±1 (see Appendix [BJ) . In the quark sector (i.e., in the 
fundamental representation), and for a generic configuration of Aq(x), no eigenvalue of will be —1 in the bulk and 
so such divergence can be disregarded. Unfortunately, in the gluon sector the situation is different since for any gauge 
configuration at least N — 1 eigenvalues of Cl(x) are necessarily unity. Therefore, the singular value v ^integer always 
appears when evaluating the adjoint trace in £2,3 and £3^. The infrared divergences are characteristic of massless 
theories at finite temperature [34tl35|. 

For v = 0, the infrared divergence comes solely from the static Matsubara mode, po — 0, in ipo. The corresponding 
integral over r has no natural scale and so the point of view can be taken that such divergences are automatically 
removed by dimensional regularization |53l ]. As explained in Appendix [Bj the integrals lt 2n without the static mode 
are given by the same expressions <|3.23(1 after the replacement v — > 1 + v in the first (^-function. The resulting 
prescription is then to use the formulas of £2,9 and £3^ with the replacements 

ifi(v) + i>{l - v) | P=0 ^(1 + v) + 4>{l - v) |~ =0 = -2 7E , 

+ ~ v) | P=0 - 1^(1 + v) + V{1 - v) | P=0 = -4C(3) , (3.32) 

to be made in the subspace f2 = 1 only, when taking the trace in the adjoint representation. One may worry that 
subtracting this subspace is not consistent with gauge invariance. This is not so. As will be discussed below, the 
periodicity of the effective action as a function of log(fi) is an important requirement. This property is not spoiled 
by the previous prescriptions. 

Alternatively, one can regulate the infrared divergence by including a cutoff function e~ m T in the r integral. The 
infrared finite modes are unaffected in the limit of small m. The static mode in tpo develops power-like divergences 
to be added to the result obtained through dimensional regularization. These terms are easily computed and are 3 

£ 2 ,m = ^--tr [llF 2 „ ± + 2£ 2 J , 
487T m p 



1 T 

£3,IR = 3 tr 

2407T 777, 
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+ Ei±FijEj± + EiFij^Ej (3.33) 



2 1 2 ^ 2 2-^2 

Even though this is a gluonic term, the result has been expressed in the fundamental representation, which is often 
preferable. (Unfortunately this is not so easily done for the other gluonic contributions, for a general SU(JV) group, 
due to the presence of the Polyakov loop in the formulas.) In these expressions we have used the notation F^ v \\ to 
denote the pieces of F^ which commute with Q, and F^± for the remainder. Specifically, in the gauge in which f2 is 
diagonal, F^ is the diagonal part of F^. As shown in Appendix [D] only terms involving at least one perpendicular 
component may be infrared divergent, and this is verified in l|3.3H|l . 

We have used here the MS scheme in dimensional regularization. Alternatively one can use Pauli-Villars regular- 
ization which amounts to inserting a regulating factor (1 — e~ rM ) in the r integration |33l |. All convergent integrals 



d Note that IT-^- a l so contains ipo and so is also afflicted by the divergence. This implies that introducing e m T is not equivalent to 

a regularization of the digamma function (and its derivatives) in the final formulas, since simple scaling relations of the type 13.151 or 
IB6I no longer hold. 
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(including 1^-) are unchanged in the limit of large M, whereas 

7 + PV = 21og(M/ M ) + 21og(^/3/47r) - ${y) - ^(1 - v) + O^r 1 ) , < v < 1 , 

7o/ V = 21og(M/ M ) + 21og(^/3/47r) - V(± + j/) - ^(i - v) + 0{M- 1 ) , -~ < v < \ . (3.34) 



(Note that these formulas do not actually depend on the scale fi.) The Pauli-Villars renormalized result is obtained by 
combining log(JVf In ) with the bare coupling constant in the tree level Lagrangian to yield the renormalized coupling 
constant g PV (/z). If, as usual, the A# parameter in the scheme R is defined as the scale fi — A# for which 1/gfjGu) 
vanishes, it is found that the Pauli-Villars and MS schemes give identical renormalized results, at one-loop, when 

MASy/Ajg) = n _2 Nf/N - ( 3 - 35 ) 

The difference between both scales comes from the yj- in (|3.26|) . which is due to the — 2e extra gluon polarizations in 
the dimensional regularization scheme |54| . 



D. Results for SU(2) and SU(3) 

We can particularize our formulas for SU(2) by working out the color traces explicitly. We use the antihermitian 
su(2) basis a/2i, so 



An 



-jS-F^, etc. 



(3.36) 



It is convenient to choose the "Polyakov gauge", in which Aq is time independent and diagonal [4||. In SU(2), 
Aq — —^iaztfi. In this case the eigenvalues of the Polyakov loop in the fundamental representation are exp(±i/30/2), 
and in the adjoint representation are exp(±i/3</>) and 1. Full results for £0,2,3(2;) in both sectors are given in Appendix 
mi Here we quote the results for £2(2) from the gluon and quark loops: 



487T 2 



( 2 log (4^) - ^ +v) - ^ Ei 



(21og( 



4ttT 



(3.37) 



with v = (/30/47T + 1/2) (mod 1) — 1/2 and = \tij],Fj k is the magnetic field. 



£2,3(2:) 



Here v — /3(/>/2tt (mod 1), and 



11 

"48^2 



21og(/i/47rT) 



1 

IT 



i± 



21og( M /47rT) 



11 



nT 



21ogW47rT) 



Ei — Eq + Ei± , Bi — Bq + Bi± , 



(3.38) 



(3.39) 



are the decompositions of the electric and magnetic fields in the directions parallel and perpendicular to Aq. Such a 
decomposition is gauge invariant. 

The quark and gluon sector contributions are periodic in <fi with periods 47rT and 2itT respectively. This periodicity 
in Aq of the coefficients multiplying the local operators is a consequence of gauge invariance. Indeed, after choosing the 
Polyakov gauge there is still freedom to make further non stationary gauge transformations within this gauge. Such 
transformations (named discrete transformations in |45j ^l are of the form U(xq) = exp(iroA), where A is a constant 
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diagonal matrix. Its eigenvalues Xj,j = l,...,N (we consider a general SU(7V) group in this discussion) are quantized 
by the requirement of periodicity in xq. For quarks, U(xq) must be strictly periodic and hence Xj = 2irinj/ (3 1 rij G Z 
(the integers rij are ^-independent by continuity). Since under a discrete transformation Aq(x) — > Aq(x) + A, the 
eigenvalues of log(fi) / (2ni) change as Vj — > i/j — rij. In SU(2) this implies that the effective action in the quark 
sector must be periodic in <f> with period AmT . In the gluon sector, periodicity of A^{x) in xq only requires that 
U(x + J) = e 2wik / N U{x ), k = 1, . . . ,N and there is an additional symmetry associated to the center of the gauge 
group 0)MIE3- That is, Xj — 2-ni(rij + k/N)/(3 in the absence of quarks (note that k is both ^-independent and 
j-independent). The eigenvalues of log(il)/(27ri) change as i/jt :— Vj — Uji — rij + rig and the effective action in 
the gluon sector must be invariant under such replacement. In SU(2) it corresponds to periodicity in 4> with period 
lixT . From this discussion it follows that an expansion in powers of log(17) breaks gauge invariance under discrete 
gauge transformations. The local operators E %, B %, Ef ± , and Bfj_ are directly gauge invariant. 

We can compare these results with those in [33 . That work goes beyond ours in that we compute the lowest terms 
in an expansion in Do whereas in |33| all orders in Aq are retained in the electric sector. On the other hand, unlike 
|33|. we treat groups other than SU(2), our gauge field configurations are not stationary, we consider higher order 
terms in the spatial covariant derivatives and we include the quark sector. 

Let us restrict ourselves to stationary gauge configurations and the gluon sector in SU(2), as in |3^|. In a notation 
close to that in |33j . the terms of the effective Lagrangian which are quadratic in F^u, but of any order in Aq, are of 
the form 

- M<t>)E\ - fi(</>)Ef ± - hz{<t>)B\ - fti($3?x • (3.40) 

To obtain these SU(2) group structure functions in our expansion we would need to retain terms with two or four 
spatial indices but any number of commutators [Aq, ]. Nevertheless, in the parallel space our calculation is complete 
since all terms of the form (DqF^u, n > 1, vanish identically in the stationary case. This implies that fs{4>) and 
h 3 ((j)) do not get any further contribution beyond those in £2,5(2), and indeed, after passing to the Pauli-Villars 
scheme with Apy = e 1 / 22 A S jg, one verifies that / 3 and h 3 of |33| are reproduced. 4 fi is not reproduced to mass 
dimension 6, but h\ is reproduced when we retain mass dimension 4 terms only, since in the magnetic sector the 
calculation in introduces ad hoc simplifying approximations which in practice are equivalent to using £,2, g {x). 

An important point is that of the periodicity of the structure functions, also emphasized in j^. In our calculation, 
the coefficients of the local operators will always be periodic in <p due to gauge invariance. Yet, this does not imply 
that the structure functions themselves should be periodic. The ones in the parallel sector, which coincide to all 
orders with the coefficients in £2,9(2), will certainly be periodic, but fi and hi will not be periodic in cj>. For instance, 
hi receives a contribution from £3^(2;) of the form /(^B^ (see Appendix O) . The function /(</>) is periodic and 
so this contribution is fully gauge invariant. However, the operator /(</>)i?Q i has still to be brought to the standard 
form in the (|3.4()|l . Using i?oi = Aq x Bj 1 , it follows that hi picks up a gauge invariant but non periodic contribution 
4> 2 f(4>). (At this point we disagree with [33j which notes that fi needs not be periodic but requires periodicity of hi.) 
We also note that in our calculation, /1 and hi are both infrared divergent, whereas in the calculation of [33 only 
hi is divergent. This should indicate that a resummation to all orders in Dq of our expansion may remove spurious 
infrared divergences. 

For SU(3) we present explicit results for the effective Lagrangian up to mass dimension 4 included. We use the 
convention 

A = -~X S A \ = -~X-A , = ~X S F° V , etc (3.41) 

where A s , s = 1, . . . , 8, are the Gell-Mann matrices. In the Polyakov gauge 

A = -iy h - t^-Xsfo ■ (3.42) 
The effective Lagrangian from the quark sector can be expressed in terms of the quantities 

"1 = J^^s + <h)> v 2 = —(-03 + <Ps) , f 3 = - ^f^8 (3-43) 



4 After correcting an inaccuracy in the Pauli-Villars treatment of {3 in V\A . The correct PV result contains 2 log(Apv/47rT) — 2/11 instead 
of just 21og(A PV /47rT). 
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as 



and 



4>,, = -^^(-| + (l-4^) 2 + (l-4l^) 2 + (l-4^) 2 ) (8.44, 



-2,9 



24tt 2 



(/" (^i) + /"N) ((^) 2 + (Flf + (Fp 2 ) 



12(4tt) 2 v v ' J ^» v > v ""' 



12(4tt) 2 v v l > J v J/; vv v 

-12^5 (/-(«*) + rw) ((^) a + (O 2 ) 

AT. 

' (f-(vi)-f-(v 2 ))F3 u F*„, (3.45) 



6V3(4tt) 2 



where we have defined 



/"W = + I7=(i/+i)(modl)-±. (3.46) 



In the gluon sector, we introduce the invariants 



"12 = ^f^ 3 ' ^31 = -^f(^3 + 30 8 ) , ^23 = ^;(-03 + 30 8 ) , (3-47) 



in terms of which the effective Lagrangian is 



C Jx) = |tT 2 T 4 (~ + p2 a(1 _ £ i2 )2 + P 2 i(1 _ £^2 + vi . , _ ;}2;i) 2j ( , JSj 
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(V(0) + f + (u 12 ) + \f + M + \f + M^j ((Flu) 2 + {F%?) 
(V(0) + \f + (vi2) + />3i) + i/ + (^ 23 )) ((Fp 2 + (F* v f) 

(/ + (o) + lf + M + \f + M + /> 23 )) ((ft) 2 + (Flf) 



' 1 1 f 2f+( Vl2 ) + V^l) + V(^23) ) (C) 2 



(4tt) 2 12 V"'' V "" 7 ' 2 J v "° 17 ' 2 
,1 (/ + (^3i) + / + W)(^) 



with 



(4tt) 2 u 

WiT! (/>3l) - /+( " 23)) ' (3 ' 49) 



= ^(P) + V(l - (v^Z), 9 = ^ (modi), 
/+(0) = -2 1E . (3.50) 
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Finally, the renormalized tree level is 

Arcc(x) = 4i^)^- (3 ' 51) 

In the stationary case, the most general structure compatible with SU(3) symmetry, constructed with two E^s and 
any number of Aq's, contains six structure functions (see Appendix lD|l 

/ 12 (&, &) {{Elf + {Elf) + fafa, <M {{E?f + {Elf) + f 67 (h^s) {{E%f + {Ejf) 

+f 33 {h,M {Eff + fss(<h, 0s) {Elf + / 38 (03>8) {EfEf) . (3.52) 

(And similarly for BiBi, etc.) Our results for £ 2 are of this form. Our expressions corresponding to /33, f$$ and /3s 
are already correct to all orders in Aq, since all Dq operators cancel in the directions 3 and 8 of the adjoint space. More 
generally for any SU(iV) and any structure function, A® decomposes F^ v into a parallel component (that commutes 
with Aq) and a perpendicular component (fully off-diagonal in the gauge in which Aq is diagonal). The structure 
functions not involving perpendicular components depend periodically on Aq and can be computed exactly using the 
appropriate finite order of our expansion (that is, the lowest order at which the corresponding local operator appears 
in£). 

In Appendix IdI we give further details on the calculation for SU(3) and SU(JV). 

IV. THE DIMENSIONALITY REDUCED EFFECTIVE THEORY 

As is well-known, in the high temperature limit non stationary fluctuations become heavy and are therefore sup- 
jCD to behave as an effective three-dimensional theory for the stationary configurations only 
Our previous calculation of the effective action was obtained by separating background 
from fluctuation and integrating the latter to one-loop. Clearly, we can adapt that procedure to obtain the action 
of the dimensionally reduced effective theory, to be denoted £'{x), by i) using stationary backgrounds and ii) taking 
purely non-stationary fluctuations only, that is, removing the static Matsubara mode in all frequency summations. 
In addition, there is a further factor (3 in C'(x) from the time integration. Note that C'(x) is not the effective action 
(or Lagrangian) of the dimensionally reduced theory but its true action (within the one-loop approximation) , in the 
sense that functional integration over the stationary configurations with £ '(x) yields the partition function. Besides 
taking A^ stationary, we will assume that Aq is small (in particular \u\ < 1), which is correct in the high temperature 
regime. We will come back to this point later. 

The static Matsubara mode is not present in the quark sector, so for that sector we simply find C' q (x) — (3C q {x). 
Likewise, the removal of the static mode is irrelevant in the ultraviolet region, hence C' tTee (x) = pC ty:ee (x) for the 
renormalized tree level. 

As discussed in Appendix 151 the removal of the static mode in the one- loop gluon sector (and for < 1) corresponds 
to replacing £(l + 2£ + 2e, v ) -> £(l + 2£ + 2e, l + i/) in lpQ3)l . For the effective potential this means B A (p) -> B 4 (l+V) 
in H3.24fl . and so (dropping an Ao-independent term) 

C'oJx) = ^T 3 tr[d 2 (l+D 2 )] , £ = log(n)/(27n), -1<P<1. (4.1) 

The analogous replacement in the mass dimension four and six terms gives (using the identity if)(l-\-p)-\-tp(l — u) = 
i){v) + i/)(-v) ) 

4,» = 7Z^ rr [§ (siog'W^T) + i - V(z?) - El - ±Ef\ , (4.2) 



£ 3, s (*) - 2 (47T) 4 T3 tr 
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-^■F^FvxFx^ + -F XfiV - 7^E^ V + -F 0fiV - j^E u + —j 



In these expressions Do stands for [Aq, }. Note that, having removed the static mode, C'(x) is free from infrared 
divergences. 
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At high temperature the effective potential suppresses configurations with f2(a;) far from unity, so by means of 
a suitable gauge transformation we can assume that A (x) is small 5 . In the absence of quarks, the situation is 
similar although in this case f2(cc) lies near a center of the group element; the center symmetry is spontaneously 
broken signaling the deconfining phase 0, IllL l55j . After a suitable generalized (many- valued) gauge transformation 
the configuration can be brought to the small Aq(x) region. It can be noted that only when Aq is small (\u\ < 1) the 
non static fluctuations are the heavy ones. If we were to choose the gauge so that v is near some other integer value 
n, the light mode would be the n-th Matsubara mode and integrating out this light mode would yield a non-local 
(and so non useful) action for the effective theory. 

Because A is small, it is standard to expand the £'(x) in powers of A Q , using the relation v = —A /{2iriT) either 
in the fundamental or the adjoint representations. We expand up to and including terms of dimension six, where now 
A coming from Q counts as dimension one. Note that this new counting is free from any ambiguity (although it is 
in conflict with explicit gauge invariance). 

The effective potential is already a polynomial in Aq. From (|3.11() and ^4.1(1 . we obtain 

We have introduced the short-hand notation (X) := tr(X) (trace in the fundamental representation) and used the 
SU(JV) identity (|D6|) . This result agrees with |39ll42l| (there written in the adjoint representation). 
In particular for SU(2) and SU(3), using the identity <D7|) valid for those groups, we find 

C' {x)=-^j + ^T{A 2 ) + ^ f (6 + N-N f ){A 2 ) 2 , N = 2, 3, (4.5) 

which reproduces the result quoted in ^} an d EJ for JV = 3. We note that consistency requires to include up to 
two- loop contributions in the effective potential [57]. 

The terms of dimension four with derivatives come from £ 2 (x), given essentially in (|3.3U|) (with -0(1 — v) —> ip(—v) 
and an extra factor /3), and setting V and v to zero. The result can be written as (the subindex 4 indicates operators 
of dimension four, all gluon fields count as mass dimension one) 

(Once again in the fundamental representation.) For the (chromo)electric and magnetic effective couplings we find 
i|W = m ~ WogW^T) + 1E) + (N + 8Nf (log 2 - \ 

' ' 2/3 (log(/i/47rT)+ 7B ) + ^ w (-iV + 8Ar / log2) (4.7) 



->M 



(T) g2( M ) ^ — > 3(4tt) 



It is possible to rescale A, and Aq (with different renormalization factors) so that C,^ (x) looks like the zero temper- 
ature renormalized tree level (|3.28|) 38, However, we will work with the original variables. 

The result for gj^(T) coincides with for A^ = 3. It also agrees with (setting Nf = 0) assuming a suitable 
Af-dependent factor between the scales A there and \x here. The scale independent ratio 



-%(T) _ 2g 2 Qu) 
- 2 r (T) 3(4^ 2 



= l-^j^(N-N f ) + 0(g^, (4.8) 



found here differs from that reference. On the other hand, in analogy with 

1 



2 M 



2/3 log(/i/AMs)> ( 4 - 9 ) 



5 To bring Aq to the \u\ < 1 basin it will be necessary to use a discrete gauge transformation, as described in the paragraph after 13.391 . 
Because such transformations are global (a;-independcnt) this will be only possible if the original Aq(x) lies in the same basin (i.e., near 
the same i nteg er u) for all x. We assume this, since otherwise Q(x) would be far from unity in the crossover region, thereby increasing 
the energy l- r >(i . 
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magnetic and electric thermal A parameters can be introduced [4 

1 



2/3 log(T/A^ M ), 



which set the scale of high temperatures for both coupling constants. For the magnetic sector we find 

N - 8N f log 2 



logCAjf/Aj^s) = 7 B - log(47r) 



22N - 4JV 



(4.10) 



(4.11) 



in agreement with [43l |. 

Next, we consider terms of dimension six. They come from C' 2 {x) expanding the digamma functions to second 
order in v and from C 3 (x) to zeroth order. From the quark sector we obtain 
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/3 3 



c (6u( x ) = i^ 3 ) (fc* ^ \ F ^ F ^ F ^ + 6i w- + 2 F °^ + 30j4 o F A " 3 ^ + QE i F v E iJ > ( 4 - 12 ) 

where we have made use of the identity F x — 2F 2 puV — AF lLV F v \F\ yi , valid inside the functional trace 0]- For gluons 
we have instead 



P 3 



57 



C UJ X ) = _ 45^ 3 )(CT tr [ F *» F *x F >v + ~2 F l^ + 27F o% + 165Agi^ v - 3£g + GE^E 
Using l|D5(l and (|D6fl . this gives for the full result 



(4.13) 



2 C(3) 
15 (4tt) 4 T 3 



(|jV - ^N f ) (F^F^F^) + (19N - 28N f )(F^) + (18N ~ 21N f )(F^) 



+(1WN - U0N f )(AtF*) - (2N - UN f )(tf*) + (4iV - 28N f )(E i F ij E J ) 



+110(A 2 ){Fl) + 220(A F tlu ) 2 



(4.14) 



For SU(2) and SU(3) the term with (AqF% v ) can be eliminated by using the identity (|D7|I . In addition, in SU(2) 
the term with (Fq„ u ) can also be removed using (|D8|> . This produces 



£(6) ( x ) 



2 C(3) 
15 (4tt) 4 T 3 
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(3 - 77V/) -F MV F vX F Xll - -Fq - 2Ei + AE^E, + (57 - 28N f ){F 2 ) 



for N = 3 , 



165 - —N f {(Al)(F^) + 2(A F^) 2 ) 



(4.15) 



^(6) ( X ) — 



4 C(3) 
15 (4tt) 4 T 3 



(2 - 77V/) ( -F^F^F^ - E* + 2E i F ij E j ) + (19 - UN f )(F 2 ^) 



-(74 - UN f ){A 2 )(F* ) + (146 - 217V / )(A F^) 2 



for N = 2. 



(4.16) 



jC'^ 6 ^(x) has been computed previously in |42j for the gluon sector and arbitrary number of colors. Our result agrees 
with that calculation (and disagrees with 0]). The dimension six Lagrangian in the quark sector has been computed 
in 0] f° r SU(3), in the absence of chromomagnetic field (A; = 0) and neglecting terms with more than two spatial 
derivatives (i.e., neglecting EfA. Our result reproduces that calculation in that limit as well. 



V. CONCLUSIONS 



In the present work we have developed in full detail the heat kernel expansion at finite temperature introduced 
in [53. We have paid special attention to the role played by the untraced Polyakov loop or thermal Wilson line in 
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maintaining manifest gauge invariance. This is a highly non trivial problem since preserving gauge invariance at finite 
temperature requires infinite orders in perturbation theory. The conflict between finite order perturbation theory and 
finite temperature gauge invariance has been previously illustrated e.g. in the radiatively induced Chern-Simons action 
of (2 + f )-dimensional fermionic theories [isj. In the case where the heat bath is chosen to be at rest the Polyakov loop 
is generated by the imaginary time component of the gauge field and can be regarded as a non-Abelian generalization 
of the well-known chemical potential. Actually, we have provided arguments supporting this interpretation; if the 
Polyakov loop was absent or represented in perturbation theory the particle number could not be fixed, as one expects 
from standard thermodynamics requirements. The new ingredient of our technique is that a certain combination of 
the Polyakov loop and the temperature has to be treated as an independent variable, in order to guarantee manifest 
gauge invariance. This can be done without fixing the gauge. 

An immediate application of our method can be found in QCD at finite temperature in the region of phcnomcno- 
logical interest corresponding to the quark-gluon plasma phase. In fact, the heat kernel expansion corresponds in 
this case to a high temperature derivative expansion organized in a very efficient way. In the case of QCD the finite 
temperature heat kernel expansion can be applied to compute the one-loop effective action stemming from the fermion 
determinant and from the bosonic determinant corresponding to gluonic fluctuations around a given background field. 
As a result we have been able to reproduce previous partial calculations and to extend them up to terms of order 
T~ 2 including the Polyakov loop effects, for a general gauge group SU(JV). As a by product we have computed the 
action of the dimensionally reduced effective theory to the same order. Further we have studied the emerging group 
structures in the case of two and three colors. 
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APPENDIX A 

Let us establish the commutation rules (|2.36(l . It is sufficient to consider the case [X,f] since D^f is a particular 
case. Because / is a function of il, it is also function of Dq through the relationship £1 = er^ Da . In fact, it is better 
to prove the relation for a general /(-Do) ( n ot necessarily periodic in its argument). No special property of Dq is 
required, so the statement is that, for any two operators X and Y and for any function /, 

[X, f(Y)} = -f'(Y)[Y, X] + \f"{Y)[Y, [Y, X}} + ■ ■ ■ 

= Y / [ —rf (n) (Y)U Y (X), D Y :=[Y,}. (Al) 

n=l 

It is sufficient to prove this identity for functions of the type f(Y) = e~ XY , where A is a c- number, since the general 
case is then obtained through Fourier decomposition. The r.h.s. of IjAlfl is 

£ ^je- XY D Y (X) = e- XY (e XD - - 1)X = e~ XY (e XY X e~ XY — X) = [X, e~ XY ] , (A2) 



n—l 



which coincides with the l.h.s. of IjAlfl . We have used the well-known identity e Dv (X) = e Y X 



e~ Y . 



APPENDIX B 

The basic integrals are 

/>oo 

Z±(i/,a):=/ drT Q -V^(e 2 ™), i/.aeR, n = 0,l,2,... (Bl) 
Jo 

where the functions (p n are defined in (|2.15(l and ± refers to the bosonic and fermionic versions, respectively. For the 
bosonic version. 



# (*, a) = ^f. (*f) n $> - v? f dTT ^ y , e -m 2 ^ t „f z . (B2) 
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We have excluded the case e 27TW = 1 which is discussed below. Integration over r gives 

n r(a + | + |) (0\ 2a y, {k-vY i 



Defining j/ = fc + < z/ < 1, the sum over k can be split into the sum for k < fco and another for k > k . In terms 
of the generalized Riemann ^-function [5lJ this gives 

r( | } 2 (^j [(~«) n C(2a + l,P) + z"C(2«+l,l-P)] , 0<?<1, ^ = fc + £, fc GZ.(B4) 
For the fermionic version, using ip~ (w) = (/?+(— w) (and so v — > t/ + i), one obtains 

4" «) = ^^(l^^ (^) " [(~ l )"C(2a + 1, 1 + V) + f\(2a + 1, ± - F)] , ~ ^ < ^ < ^ . (B5) 

Note that 

4,(1/, a) - (-lr ^+^Jh t^a). (B6) 

The formulas are consistent with periodicity and parity 

I±(v,a) = fi{v + l,a) = (-l) n lt(-v,a) . (B7) 

As discussed in IIVI the dimensionally reduced effective theory for the stationary configurations, requires to remove 
the static mode from the summation over Matsubara frequencies in the bosonic integrals. This prescription breaks 
periodicity in v but this is not relevant for the effective theory, since it only describes the small A§ (or v) region. (A 
prescription that preserves periodicity would be to remove the frequency k = k when v < ^ and k = k + 1 when 
v > |.) The result for the \v\ < 1 is 

«) = r( i } 2 [i) [H)"C(2a + 1, 1 + v) + *"C(2« + 1, 1 - 1/)] , -K 1/ < 1 . (B8) 

A related issue is that of the infrared divergences for integer v. Due to periodicity, we can restrict the discussion 
to the case v = 0. For n ^ 0, the static Matsubara mode does not contribute to I+(i/,a), and so there is no infrared 
divergence in this case. On the other hand, in Iq(v, a), the static mode is either infrared or ultraviolet divergent. In 
dimensional regularization such an integral [y = k = n = in ljB2|) ~) is defined as zero since it has no natural scale 
|53|. So for all n the result is equivalent to removing the static mode 

/>, «) = 4 + (0, a) = { (-D" /2 2-- 1/2 r(a + i + 03/2-) 2 " C(2a + 1) , even n for y g z (Bg) 

Alternatively one can regulate the infrared divergence by adding a cutoff function e~ m T (m — > 0) in the r integral. 
This amounts to adding a contribution v / 47rT(a + l/2)/(/3m 2a+1 ) in 2q~(2/, a) for integer v. 

APPENDIX C 

In this appendix we present results for SU(2) in both sectors, including all terms of mass dimension 6. All results 
are given in the MS scheme. In these formulas we have allowed for an explicit infrared cut off m, as commented 
at the end of Appendix [Bj The results with strict dimensional regularization are recovered by removing all infrared 
divergent terms from the formulas. The conventions are those of subsection IIII Dl 

Aree(^) = 73 l 7 - T ^, (CI) 



4g 2 M 
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£o, fl (*) = ^ (-i + 4P 2 (l-9) 2 ) , (C2) 
CM = -g^(n + 2Io «(i^)-^-^ 1 -^)^ll 

+ ^~M^)^ (C3) 



1 / 1 ^ 2 



(<//'(£) 9)) F 2 



1 / 1 \ 2 / o ( nT\ 3 



96tt 2 V4ttT, 
1 / 1 ^ 2 



80tt 2 V47rT 



1 ^ 2 



(V"(^)+V"(1-^))4V|I 



; V/ ^^V + 4c(3)-r"(^)- ( '"(:i-^))^ 



1 / 1 \ 2 /ttT x 3 



IOtt 2 \4itT m 0l± 



2 



+ 24^(^t) «(^)+<a-))4„ 



1 / 1 x2 



240vr 2 V 4ttT 



+ ^"(1 - P)) Eijfc^i x Ej) ■ B k 



+ 24^ (i) 2 ( 8 ) 3 - 4C(3) ^" {d) - 9) ) e «*^ X ^ ■ ' (C4) 

Cojx) = 2 r 2 T±N f (1 1(1 - 477 2 ) 2 ) , (C5) 



x [—{F^ x F^) • + -F^ v - F^ u - 2e ijk (Ei x Ej) ■ B k + 3F^ U - 2E? i j . (C7) 
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It can be noted that the quark terms do not distinguish between parallel and perpendicular components. This is 
due to the fact that in SU(2) an even function of 77 (or any other element of su(2)) in the fundamental representation 
is necessarily a c-number. Since the ip n functions involved to mass dimension 6 are all even, the V dependence gets 
out of the trace in (|3.14|) and i|3.16fl and Ao is no longer a privileged direction in color space. This mechanism does 
not act in the adjoint representation, i.e., in the gluon sector, nor for other SU(iV) groups (cf. (|3.45(l ). 

The infrared divergence is tied to v integer, so it does not exist for fermions, and also cancels in all gluon terms 
involving only parallel components. 



APPENDIX D 



In SU(iV) the gauge can be chosen so that Aq is diagonal. This form is unique (up to permutation of eigenvalues) 
and produces N — 1 quantities invariant under SU(^V) (^3, for SU(3)). If X represents or any other element 
of su(iV) (X* = —X, tr(X) — 0) with N 2 — 1 independent components, we can use the remaining gauge freedom (the 
N — 1 gauge transformations which leave Ao diagonal) to fix N — 1 of these components. This adds (TV 2 — 1) — (N — 1) 
new invariants involving X (and Ao). Of these, N — 1 are linear in X (the diagonal components of X), N(N — l)/2 
are quadratic and (N — 1)(N — 2)/2 are cubic. For instance, in SU(3), under a diagonal gauge transformation 

b \ / x e< a -&a e^ 2a+ ^b\ 

c ^ _ e -*(a-/J) * y e ^+2(i) c ( D1 ) 

X-yJ \^ e -m*+f}) h * _ e -*(a+2/3) c * - X -y j 

the invariants are x, y, aa* , bb* , cc* and ab*c (the last one is complex but its modulus is not independent). For X = E{ 
this gives the six structure functions in (|3.52|) . Each further vector Y £ su(N) produces new N 2 — 1 invariants. 

For computing the traces in the adjoint representation one possibility is to use the adjoint basis (T s ) rt — f rs t, such 
that to F^u — F^ v t s (t s — \ s /2i) in the fundamental representation it corresponds F^ = F^ V T S in the adjoint one. 
We have also used an alternative approach, as follows. The elements of su(iV), such as the gluon quantum fluctuation 
are N x N matrices, (a^ad- From the action F flI/ (ax) = [F^, a\], it follows 

(F^)aa,bb = ( F ^ S ab- S ab(F^) U , a,b,a,b = 1, . . . , N . (D2) 

In matrix notation this can be written as F^ — F^ u <£> 1 — 1 ® F^ v — F^ u ® 1 + 1 ® F* u , or even, in shorter form, 

F^u = Ffi V — F^ v = F^ v + F^ u (D3) 

understanding that F^ v or F* always refer to the dotted space. Similarly, = A^ — A T ^ = A^ + A*. Since 
dotted and undotted operators commute, it follows that the f2 = fig)fl* = fi<8) il~ 1T for the Polyakov loop. In 
the Polyakov gauge (Aq stationary and diagonal) Q is diagonal (Q) a b = u a 8 a b and VL is also diagonal in that basis, 

(tyaa,bb = Uaa5ab5abi with W «i = W aW['. 

From the point of view of the gauge group, the computation of the trace in the adjoint space involves only four 
different structures appearing in &J g ,&^ g , and bj g . These are 

trCf(fi)) = ^'/( Wm)i 
aa 

tr(/(0)l^) = /W [(Fpaa + (Fpaa - 2(F ltv ) aa (F lo ,) ii ] , 

ad 

t r (i '{tyF^vFvxFx^)) = f(w a a) [(FfivFvxFx^aa + (F^F^xFxfi) dd ~ {F^aaiF^xFx^da ~ (F^F l/ x)aa(Fx t i)dd\ , 

ad 

te{f{&)EiFijEj) = f(^ad) [(EiFijEj) aa + (EiFijEj)ad — (Ei) aa ([Fij , Ej]) aa — ([Ei, Fij]) aa (Ej) aa 

ad 

— (EiEj) aa (Fij)aa — {Fij) aa {EiEj)aa\ . (D4) 

(52'ad m t ne nrs t equation indicates that one of the N modes with a = a should not be included. This removes the 
singlet mode present in U(N) but not SU(iV). The singlet mode does not contribute in the other formulas.) Often, 
f(u>) = /(lu^ 1 ) (i.e. f(uj a a) is symmetric in a, a), but this property has been not used here. It can be observed that 



Cx a 
-a* y 
-b* -c* - 
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the contributions a — a, which correspond to 51 = 1 and are afflicted by infrared divergences, cancel in the subspace 
parallel (i.e. for F^ u diagonal in the Polyakov gauge). 

Useful SU(iV) identities are (( ) stands for trace in the fundamental representation) 

tr (X 2 ) = 2N(X 2 ) , X G su(JV) , (D5) 



tr(X 2 Y 2 ) = 2N(X 2 Y 2 ) + 2(X 2 ) (Y 2 ) + A(XY) 2 , X, Y G su(iV) 



(X 2 Y 2 ) = -\{[X,Y] 2 ) + \{X 2 ){Y 2 ) + 1 -{XY) 2 X,Ye su(3) , 



([X, Y] 2 ) = -2(X 2 ) (Y 2 ) + 2(XY) 2 X,Ye su(2) . 



(D6) 
(D7) 
(D8) 
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